Role of the superconducting gap opening on vertex corrections 
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The nonadiabatic effects due to the breakdown of Migdal's theorem in high-T c superconduc- 
tors are strongly affected by the opening of the superconducting gap. Here we report how the 
momentum-frequency dependence of vertex corrections is modified by the gap. A general effect is 
that the positive region of the vertex corrections is increased leading to an enhancement of their 
relevance. This has a number of physical consequences that we discuss in details in relation to 
specific experiments. 
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I. INTRODUCTION 

A common peculiar characteristic of many unconven- 
tional high-Tc superconductors (cuprates, A^Cqq com- 
pounds,...) is the narrowness of the electronic bands 
crossing the Fermi level, leading to Fermi energies tp of 
the same order of magnitude of the phonon energies u> p h 
{ujph/tF^- 0.2 - 0.3) .EJ In such a situation, the adiabatig 
assumption {oj v h/£f -C 1) on which Migdal's theoremel 
rests cannot be used anymore to justify the omission of 
vertex corrections in a diagrammatic theory. Motivated 
by this observation, a renewed interest has recently arose 
about the possible effects due to inclusion of the vertex 
corrections in the electron-phonon (or any kind of bosonic 
mediator) interaction. This task is even more relevant 
for the purely electronic interaction models (Hubbard, of 
t — J), whem no justification exists at all to neglect such 
corrections .13 

In this perspective, in previous papers, we have per- 
formed an intensive study of the vertex corrections in 
the normal state, and, afterwards, we have generalized 
the conventional Migdal-Eliashberg theory in order to in- 
clude the first nonadiabatic-<cp|rrections due to the break- 
down of Migdal's theoremcO We have investigated the 
effect of the Migdal corrections on different quantities re- 
lated to both the superconducting and the normal state 
properties, as for instance the superconducting critical 
temperature, or the isotope effect on the effective electron 
mass.Q It was shown that an important feature of these 
corrections is the non trivial structure in momenta and 
frequencies space, so that the actual relevance of their 
inclusion can, or cannot, affect in a significant way the 
properties of the system depending on the effective pa- 
rameters of the systems. In particular, the critical tem- 
perature T c is found to be sensitive to the momentum 



transfer of the electron-phonon interaction, and it would 
be interesting to evaluate the effect of the vertex correc- 
tions to the gap-T c ratio. Moreover, a negative isotope 
effect is predicted on the effective electron mass to* in 
the normal state. This is in quite good agreement with 
recent experiments that show a large isotope effect on 
to* .BO However, this result is based on measurements 
performed in the superconducting state, where a nonadi- 
abatic generalization of Migdal-Eliashberg theory to in- 
clude vertex corrections is still missing, and no prediction 
on the effect of vertex corrections can be yet done .1121 

As a first step to overcome this lack, we are going to 
analyze in some detail, in the present paper, the ver- 
tex function and its momentum-frequency dependence in 
the superconducting state. Our conclusions show that, as 
soon as a superconducting gap is opened, the momentum 
and frequency structure of the vertex corrections is dras- 
tically changed, driving the system towards an effective 
strengthening of the electron-phonon coupling. 

The paper is organized as follows: in the next sec- 
tion we introduce the vertex function in the normal state 



while in Sec. Ill we study the modifications induced by 
the opening of a superconducting gap. Section IV is de- 



voted to a general discussion on the possible implications 
on physical quantities. 



II. THE VERTEX FUNCTION IN THE NORMAL 
STATE 



The lowest order correction due to the breakdown of 
Migdal's theorem in the normal state is diagrammati- 
cally depicted in Fig. 1, which defines the vertex function 
r(q,w„;k,w m ): 
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Here, e(p) and ui(k — p) are electronic and phononic 
dispersions and g(k — p) is the matrix element of the 
electron-phonon interaction, T is the temperature and 
lu u , u> m , and coi are Matsubara frequencies. 

As it was discussed in Ref. [|, the qualitative behavior 
of the momentum- frequency structure of T(q, w n ; k, w m ), 
is well caught by two characteristic limits, the static and 



the dynamical one, T s and T^, defined as follows: T s = 
r(q -> 0, uj n = 0; k, uj m ), T d = T(q = 0, w„ -> 0; k, w m ). 
In fact, it can be show that the difference between these 
two limits is in geeeral finite and it is given by the fol- 
lowing expression :til 
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An analytic calculation of the vertex function T, valid 
for small q, was also derived by considering an Holstein 
model, with Einstein phonon dispersion w(q) = lo q , a 
constant electron-phonon bare vertex g(q) = g, and a flat 
electron density of states for a band of width E. Within 
this model, at zero temperature and for a negligible ex- 
ternal frequency u> m , Equation (0) simplifies to: 



r d - r s = a, 



(3) 



where A = 2g 2 N(0)/uo is the electron-phonon coupling 
constant and N(0) is the electronic density of states at 
the Fermi level. The two limits read 



T d = A/(l + 2^/E) 
T s = -A/(l + E/2u! ) 



(4) 



where 2loq/E represents, in this model, the adiabatic pa- 
rameter of Migdal's theorem. 

The sign of the vertex function is plotted in Fig. 2 as 
function of the exchanged momentum q and frequency 
Lo n . The non-analyticity of the vertex at (q = 0, w„ = 0) 
is here clear, where the negative static limit is recovered 
on the vertical axis, and the positive dynamical limit on 
the horizontal one. This complex structure in momen- 
tum and frequency has important consequences on the 
role of the vertex corrections for different quantities. For 
instance, the renormalized electron mass to*, and so ijt 
isotope effect, is mainly related to the static limit onhv 
while for the superconducting critical temperature the 
full (q, Lo n ) dependence of the vertex function is relevant 
This difference can be reflected in a stronger dependence 
of T c , and its isotope effect, on the microscopical parame- 
ters than the isotope effect on to* , as it is experimentally 
seen .Ha 

However, we would like to stress that the above pre- 
dictions are derived in the normal state. How are these 
results affected by the onset of the superconductivity? 
Despite the difficulty to build a generalization of Migdal- 
Eliashberg theory including the nonadiabatic corrections 
to the Migdal's theorem in the superconducting state, 



some preliminary answers at this question can be ob- 
tained by analyzing the modifications on the vertex func- 
tion induced by the superconducting order parameter. 



III. THE VERTEX FUNCTION IN THE 
SUPERCONDUCTING STATE 

The existence of a long-range superconducting order 
parameter modifies the vertex correction in two ways 
(Fig. 3). On one hand, the single particle propagators 
involved in the diagram (Fig. 3a) are affected by the 
opening of the superconducting gap. On the other hand, 
a new diagram appears, involving two anomalous Green's 
functions (Fig. 3b). We label them T on and T°f f , respec- 
tively, reminding that they are related to on-diagonal and 
off-diagonal propagators. In order to preserve the anal- 
ogy with the previous calculations at T > T c (Fig. 1), r°" 
and T°ff are obtained by using Green's functions renor- 
malized only by anomalous self-energies, and the normal 
state renormalization is not taken in account. Within 
this approximation, the normal and anomalous propaga- 
tors are given by: 



G(k, iwr, 
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(5) 
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where E(k) — y/e(k) 2 + A 2 is the BCS superconducting 
excitation spectrum and A is the order parameter. In this 
paper we consider only isotropic s-wave symmetry, how- 
ever most of the results will be valid also for anisotropic 
s-wave and <i-wave type of symmetries. By using the 
propagators defined above, the extension of Eq.^ in the 
superconducting state reads: 
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The first important difference of the above expression 
with respect to the vertex function in the normal state, 
Eq.(pl), is already evident by looking at the difference 



between the dynamical and static limits: 



E 



2\g(k - p)| 2 df[E(p)} uj(k - p) 2 [E(p) 2 - cj(k - p) 2 - ujI + 2iu m e(p)] 



p w(k- P ) dE( P ) [ w (k-p) 2 + <-£;(p) 2 ] 2 + [2 Wm £;(p)] 2 



(8) 



In fact, contrary to Eq. (JTJ), the right hand side of Eq.(||) 
vanishes when T — > because of the gap opening in the 
quasiparticle spectrum. One can realize such a result by 
noticing that for T <C A the integral over the Brillouin 
zone has roughly the weight factor 



df[E(p)} exp(-A/T) 
dE(p) ~ T 



(9) 



which goes to zero for T — ► 0. The vanishing of Yd — Y s 
at zero temperature is an intrinsic feature related to the 
opening of a gap in the excitation spectrum. It is easy 
to verify that such a result holds true also for a d-wave 
symmetry of the order parameter. Therefore the equal- 
ity = Y s at zero temperature characterizes the deep 
modification of the vertex correction due to the onset of 
superconducting long-range order. 

In order to compare in more details the effect of the 
superconducting gap opening on the vertex function we 
consider the same model discussed in Ref. ^, that we re- 
mind here briefly: the phonon dispersion is taken as a 
Einstein one w(q) = u>q, the electron-phonon matrix ele- 
ment constant <?(q) = g, the electronic density of states 
flat: 



N(e) = N(0) 



E E 
2 <£< 2 



where E is the total bandwidth and E/2 = Ep the Fermi 
energy. Moreover, since we are considering the vertex 
corrections at the first order in the adiabatic parameter 
cuo/Ep, the .-sum on the internal momentum p can be 
replaced by:ta 



f E/2 

V->iV(0) / cfe«...» FS 

„ J-E/2 



For the same reason, e(p + q) can be written as 
e(p + q) ~ e(p) + vp\q\ cos 9 = e + y, and the variable y 
represents the angular degrees of freedom to be averaged 
on the Fermi surface. In this context, the vertex function 
will depend on momenta just via the dimensionless vari- 
able Q = \q\/2kp, where |q| is the momentum transfer. 
Finally, since the vertex function T(q, ui n ; k, u> m ) depends 
weakly on the external frequency u m , we set cu m = 0. 
With these notations, the expression of the vertex func- 
tion r is: 
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where, as before, A = 2g 2 N(0)/uJo. The sum over u>i can 
be performed exactly by using the Poisson's formula at 
T = : — * dwj2-K. In this way, the static 
and dynamical limits can be analytically obtained. The 



derivation is quite cumbersome, but straightforward, and 
does not present particular difficulties. Namely, we ob- 
tain 
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where, for the sake of shortness, we have set m = 2ujq/E 
and 5 = A/u>q. Comparing Eq.(|ll|) with Eq.(^) we can 
see that the onset of the long-range superconducting or- 
der parameter leads to a drastic change of the structure 
of the vertex function. In particular, while the dynam- 
ical limit is affected in a smooth way by a finite A, the 
value of the static limit immediately jumps to the pos- 



itive value r s = Yd as soon as the gap is open. This 
result is not surprising, since similar modifications of the 
analytical properties of different susceptibilities (and the 
electron-phonon vertex can be described as a particular 
susceptibility) were already recovered £3113 

The global change of Y(Q,oj n ) due to the opening of 
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the superconducting gap is also more evident in Fig. 4, 
where we plot the sign of the vertex function in the Q—uj n 
space, as calculated directly by Eq. ([h]), without any 
small Q assumption. With respect to the normal state, 
marked by the dashed line in Fig. 4, the change of sign in 
the superconducting state is shifted towards larger values 
of Q, leading to an enlargement of the region of positiv- 
ity of the vertex function. Moreover, the non-analyticity 
of the vertex function in the normal state at the point 
(q = 0,cj„ = 0) is completely removed. For zero ex- 
changed frequencies u> n = 0, the vertex changes sign at 
a finite value Qa of the exchanged momentum, roughly 
given by Qa ~ A/wo- This feature can be attributed to 
the presence of a characteristic length which we identify 
with the superconducting coherence length £o- The inter- 
play of the vertex function and the coherence length will 
be discussed from the physical point of view in Sec. IV. 

In addition to the increase of the region where 
T(Q,Lj n ) > 0, the opening of the superconducting gap 
leads also to an amplification of the global magnitude 
of the vertex. This behavior is clearly shown in Fig. 
5, where we plot the average of the vertex function 
over the Brilkmin zone and over the relevant exchanged 
frequencies:! 



duj 



u 



7T UJ Z + LOq Jq 



2QT(Q,u)dQ. 



(12) 



When A = the positive and negative contributions of 
the vertex function nearly cancel each other giving a neg- 
ligible total effect also for intermediate values of ujq/Ef- 
The case A / gives instead an imbalance in favour of 
the positive contribution resulting into an enhancement 
of (r) > 0. This result is a consequence of both the en- 
largement of the positive sign region and at the same time 
an increase (suppression) of the magnitude of the posi- 
tive (negative) values of T(Q,uj n ). Despite of the deep 
modifications induced by the opening of gap, the vertex 
function still fulfills the adiabatic limit of Migdal's theo- 
rem lim Wo/BF ^ (r) = 0. 

The importance of small momentum transfer selection 
for the global sign of the vertex function in the normal 
state has already be outlined in Refs. There, an up- 
per cut-off Q c was introduced in order to select values of 
the momentum transfer Q < Q c . In this way the average 
procedure of Eq. ( |l2| ) was modified to restrict the integra- 
tion over Q according to the constraint Q < Q c . The so 
defined constrained average turns out to be very sensi- 
tive to Q c since it modifies the balance of the positive and 
negative regions of T(Q,uj n ). In particular, small values 
of Q c lead to an effective enhancement of (r). This sit- 
uation is quite altered by the opening of a finite gap. In 
fact, from Fig. 4 it can be realized that an upper cut-off 
Q c has less important effects on the global magnitude 
when A / 0. Therefore, in the superconducting state, 
the Q c dependence of (r) is more and more weakened by 
increasing the value of A. This feature is described in 
Fig. 6 where we compare (r) for A = and A = 0.5wo 
as a function of Q c - However, more than the magnitude 
of (r), the most significant quantity to look at is the 
sensitivity of the average vertex, that can be expressed 



as the derivative d(T)/dQ c (insert of Fig. 6). For small 
Qc's, that in our perspective is the relevant region for the 
high-T c materials, the sensitivity of the averaged vertex 
function (r) is much less pronounced for A ^ 0, implying 
that microscopical variations of Q c lead to smaller mod- 
ifications of (r) in the superconducting state than in the 
normal one. 

As we are going to discuss below, all these modifica- 
tions on the vertex structure induced by the gap opening 
can have important implications on physical quantities 
of the superconducting state. 



IV. DISCUSSION 

As we have seen in the above calculations, the onset of 
a long-range superconducting order parameter changes 
substantially the momentum-frequency structure of the 
vertex function. In particular, the region of positivity 
is enlarged with respect the normal state, and its rele- 
vance for small adiabatic parameter is increased. The 
main modifications can be sketched by the static and 
dynamical limits. The dynamical limit T c i is just slightly 
affected by the onset of superconductivity, and it remains 
always positive, whereas the static limit T s is drastically 
changed switching from a negative to a positive value. 
In particular, at zero temperature these two limits are 
strictly equal for A / 0. This behavior can be under- 
stood by considering that the dynamical limit, which is 
essentially independent of the many-particle properties, 
is related to the static lattice distortions, and as matter 
of fact it defines the attractive electron-phonon coupling. 
On the other hand, the static limit depends on many- 
body effects due to the Pauli exclusion principle, so that 
it is strongly determined by the fermionic.-aicture of the 
single-particle properties of the electrons .EJ Such a pic- 
ture is completely destroyed in the superconducting state 
where all the Fermi electrons are condensed in Cooper 
pairs, and Pauli principle is ineffective. As a result, no 
more difference between dynamical and static limits is 
found. Moreover, the presence of a finite value Qa of 
the momentum transfer Q for which at zero exchanged 
frequency the vertex changes sign is a consequence of the 
presence of the superconducting coherence length £o- 

Therefore, based on this result, which kind of physical 
predictions can be done? 

First of all, it is straightforward to expect that the ver- 
tex corrections beyond Migdal's theorem become more 
and more important in the superconducting state with 
the magnitude of the gap, and, since the positive part of 
the vertex function is increased, this leads to a strength- 
ening of the electron-phonon coupling. We remind that 
the magnitude of the superconducting gap A(0) is of 
the same-order of the phonon frequencies in the high-T c 
materials!!!, so that such a situation gives a huge increase 
of the positive part of the vertex function. As a first con- 
sequence, the zero temperature quantities, as for instance 
A(0), will feel a much stronger effect of vertex corrections 
with respect T c , and the ratio 2A/T C is expected to in- 
crease the BCS value 3.53. The modifications induced 
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by the superconducting gap on the vertex corrections 
are essentially ruled by ratio A(0)/cj , so that materi- 
als with high phonon energies, as the Cgo compounds! 1 ^, 
should show a lighter enhancement of 2A/T C with re- 
spect the cuprates, where the phonon energies involved 
are smaller .El This is actually recovered by experimental 
measurementsJiaO 

In addition we can argue that the superconducting 
properties are much more robust inside the supercon- 
ducting phase than the properties related to the onset 
of superconducting instabilities. In fact, as the posi- 
tive and negative parts of the vertex function are almost 
equivalent in the normal state, the critical temperature 
is strongly affected by the an implicit possible predomi- 
nance of forward small q scattering, that will select the 
positive part. Such a situation is naturally given by tak- 
ing in account the strong electronic correlation due to the 
on-site Coulomb interaction!!^, that is considered to be a 
crucial ingredient of the high-T c superconductors. The 
relevance of the forward scattering is strongly affected 
by microscopic parameters, like the adiabatic ratio and 
the doping. As a consequence, the critical temperature 
itself is expected to be strongly related to these param- 
eters. On the other hand, due to the opening of the 
superconducting gap, the vertex corrections involved in 
the zero temperature gap A(0) are essentially positive, 
and the small q selection given by the strong correlation 
is less relevant. This could explain the different behavior 
of A(0) and T c as function of the doping <5£3 
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FIG. 1. The lowest order vertex correction diagram. Solid 
lines are electronic propagators, dashed lines are phonon prop- 
agators. 
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FIG. 2. Positive and negative regions of the vertex func- 
tion in the normal state (Ref. 4). It was derived by a small 
Q expansion. The adiabatic parameter is 2ujq/E — 0.5. 
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FIG. 3. First order vertex diagram in the superconducting 
state. 
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FIG. 5. Momentum-frequency average of the vertex func- 
tion (F(Q,lo)) as function of the adiabatic parameter 2u>o/E 
for different gap magnitudes. A is set equal 1. 




FIG. 4. Comparison between positive and negative regions 
of the vertex function in the normal state (A = 0) and in the 
superconducting state (A = O.luo) for 2u>o/E — 0.5. 
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FIG. 6. Plot of (r) as function of the momentum cut-off 
Q c , for 2lj /E — 0.3 and A = 1. Insert: the corresponding 
derivative: d{T)/dQ c . 
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